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Abstract 

Particle oscillations following neutrino production processes are analysed within 
Feynman's path amplitude formulation of quantum mechanics. Consideration of 
the temporal sequence of production and detection events reveals an important 
contribution to the oscillation phase from the space-time propagator of the decaying 
source particle. Formulae are derived for spatial neutrino oscillations following pion 
decay (at rest and in flight), muon decay and nuclear /3-decay at rest, as well as for 
similar muon oscillations after pion decay at rest and in flight. In all cases studied, 
the neutrino oscillation phases found differ from the conventionally used one. A 
concise critical review is made of previous treatments of the quantum mechanics of 
neutrino and muon oscillations. 
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1 Introduction 



The quantum mechanical description of neutrino oscillations |^ has been the sub- 
ject of much discussion and debate in the recent literature. The 'standard' oscillation 
formula yielding an oscillation phaseQ, at distance L from the neutrino source, be- 
tween neutrinos, of mass nii and 1112 and momentum P, of: 

</'l2 = ^ \p'' , (1-1) 

is derived on the assumption of equal momentum at production of the two neutrino mass 
eigenstates. Other authors have proposed, instead, equal energies |^ or velocities at 
production, confirming, in both cases, the result of the standard formula. The latter refer- 
ence claims, however, that the standard expression for 0i2 should be multiplied by a factor 
of two in the case of the equal energy or equal momentum hypotheses. However, the equal 
momentum, energy or velocity assumptions are all incompatible with energy-momentum 
conservation in the neutrino production process . In two recent calculations 0, §] a co- 
variant formalism was used in which exact energy-momentum conservation was imposed. 
Using the invariant Feynman propagator to describe the space-time evolution of the 
neutrino mass eigenstates, the author of Ref. [0 also found an oscillation phase a factor of 
two larger than given by the the standard formula. However, on convoluting the transition 
amplitude with a Gaussian wavefunction descibing the spatial-temporal position of the 
source pion, the standard result was recovered. In Ref. |Q , a completely different formula 
was obtained for the neutrino oscillation phase following pion decay and it was predicted 
that only correlated oscillations of neutrinos and the recoiling decay muons could be ob- 
served. However, the author of Ref. ^ as well as others |11[] claimed that muon 
oscillations would either be completely suppressed, or essentially impossible to observe. 

The present paper calculates the probabilities of oscillation of neutrinos and muons 
produced by pions decaying both at rest and in flight, as well as the probabilities of neu- 
trino oscillation following muon decay or decay of a nucleus at rest. The calculations, 
which are fully covariant, are based on Feynman's reformulation of quantum mechan- 
ics [0 in terms of interfering amplitudes associated with classical space-time particle 
trajectories. The essential interpretational formula of this approach^, though motivated 
by the seminal paper of Dirac on the Lagrangian formulation of quantum mechanics |jl3|. 



and much developed later in the work of Feynman and other authors |Ti|, was actually 
already given by Heisenberg in 1930 [|^. The application of the path amplitude formal- 
ism to neutrino or muon oscillations is particularly staightforward, since, in the covariant 
formulation of quantum mechanics, energy and momentum are exactly conserved at all 
vertices and due to the macroscopic propagation distances of the observed neutrinos and 
muons all these particles follow essentially classical trajectories {i.e. corresponding to 
the minima of the classical action) which are rectilinear paths with constant velocities. 

^ The interference term is proportional to cos (f>i2 or sin^ . 
^Units with h/2Tr = c = 1 are used throughout. 
^Postulate 1 and Eqn.(7) of Ref. Q. 

■^Heisenberg remarked that the fundamental formula (1.2) must be distinguished from that where the 
summation over intermediate states is made at the level of probabilities, rather than amplitudes, and 
that the distinction between the two formulae is 'the centre of the whole quantum theory'. 
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The essential formula of Feynman's version of quantum mechanics, to be employed in the 
calculations presented below, is [|l^, |15[] : 



fi 



J:T.■■■J:{f\kl){kl\h)...{kn\^) 

kl k2 kn 



[1.2) 



where Pfi is the probability to observe a final state /, given an initial state z, and fcj, j = 
l,n are (unobserved) intermediate quantum states. In the applications to be described in 
this paper, Eqn(1.2) specialises to 0: 



Pe 7r+ 



Yl kfc)(^'fc,a;z)|z/fc,Xfe)(z/fc|7r+)(7r+,Xfc|7r+,Xo)(7r+,Xo|5'^) 

k=l,2 



for the case of neutrino oscillations and 



1.3) 



(e^l/^fc ' xnlfit, Xk){fMt\T^^) (tt^, Xk\7r^, xo)(vr+, xo|^^) 



fc=l,2 



:i.4) 



for the case of muon oscillations. Pe-n+ is the probability to observe the charged cur- 
rent neutrino interaction: z/gn —>■ e~p following the decay: tt^ f^'^i^fi, while Pe+7r+ is 
the probability to observe the decay fi^ — > e'^u^'U^, after the same decay process. In 
Eqns.(1.3),(1.4) \i'k),k = 1,2 are neutrino mass eigenstates while \fi^),k = 1,2 are the 
corresponding recoil muon states from pion decay. {I'kl'^'^), (/^fc Itt"*") and (e"*"!/!^) denote 
invariant decay amplitudes, (e~|z/fc) is the invariant amplitude of the charged current neu- 
trino interaction, {p, X2\p, Xi) is the invariant space-time propagator of particle p = u, fi,n 
between the space-time points xi and X2 and (tt"'", xqIS'tt) is an invariant amplitude de- 
scribing the production of the vr"*" by the source S.,, and its space-time propagation to the 
space-time point xq. 



The difference of the approach used in the present paper to previous calculations 
presented in the literature can be seen immediately on inspection of Eqns(1.3) and (1.4). 
The initial state[| is a pion at space time-point Xq, the final state an e~ or produced 
at space-time point xd. These are unique points, for any given event and do not depend 
in any way on the masses of the unobserved neutrino eigenstates propagating from Xk 
to X£) in Eqn(1.3). On the other hand the (unobserved) space-time points Xk at which 
the neutrinos and muons are produced do depend on k. Indeed, because of the different 
velocities of the propagating neutrino eigenstates, only in this case can both neutrinos 
and muons ( representing alternative classical histories of the decaying pion) both arrive 
simultanously at the unique point xd where the neutrino interaction occurs (Eqn(1.3)) or 
the muon decays (Eqn(1.4)). 

The crucial point in the above discussion is that the decaying pion, via the different 
path amplitudes in Eqns(1.3) and (1.4), interferes with itself. To modify very slightly 

^In Eqns(1.3) and (1.4) an additional summation over unobserved states, with different physical masses 
of the decay muon, is omitted for simplicity. See Eqns.(2.1) and (2.35) below. 

^ The pion production and propagation amplitude (tt"*", xqIS'^) contributes only a multiplicative con- 
stant to the transition probabilites. The initial state can then just as well be defined as 'pion at xo\ 
rather than \Sn)- This is done in the calculations presented in Section 2 below. 
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Dirac's famous statement^: 'Each pion then interferes only with itself. Interference be- 
tween two different pions never occurs'. As will be seen in the discussion in Section 
5 below, essentially all recent treatments of neutrino oscillations in the literature have 
attempted to contradict the second sentence in Dirac's statement! 

Because of the different possible decay times of the pion in the two interfering path 
amplitudes, the pion propagators (tt^, Xfc|vr+, xq) in Eqns(1.3) and (1.4) above give im- 
portant contributions to the interference phase. To the author's best knowledge, this 
effect has not been taken into account in any previously published calculation of neutrino 
oscillations. 



The results found for the oscillation phase are, for pion decays at rest: 



:i.5) 



and for pion decays in flight: 



'VI 



{ml - m^^)E^ cos 9^ 



where 



{ml — m'j^yE'^ cos 9^ 



Am = mi — mn 



:i.6) 
:i.7) 



The superscripts indicate the particles whose propagators contribute to the interference 
phase. Also -E^, Ei, and -E^ are the energies of the parent vr and the decay z/ and /i 
and 9,y, 9^ the angles between the pion and the neutrino, muon flight directions. In 
Eqns.(1.5) to (1.7) terms of order mf, and higher, are neglected, and in Eqns.(1.6) 
and (1.7) ultrarelativistic kinematics with i?7r,^ ^ fn-w^fi is assumed. Formulae for the 
oscillation phase of neutrino oscillations following muon decays or nuclear /9-decays at 
rest, calculated in a similar manner to Eqn(1.5), are given in Section 3 below. 

A brief comment is now made on the generality and the covariant nature of the cal- 
culations presented in this paper. Although the fundamental formula (1.2) is valid in 
both relativistic and non-relativistic quantum mechanics, it was developed in detail by 



Feynman [|1^, |TJ] only for the non-relativistic case. For the conditions of the calculations 
performed in the present paper (propagation of particles in free space) the invariant space- 
time propagator can either be derived (for fermions) from the Dirac equation, as originally 
done by Feynman or, more generally, from the covariant Feynman path integral for 
an arbitary massive particle, as recently done in Ref . [0 . In the latter case, the invariant 
propagator for any stable particle with pole mass m, between space-time points Xi and 
X/ in free space is given by the path integral 0]: 

{ im f^f ( dx dx \ , 1 , . 



^'Each photon then interferes only with itself. Interference between two different photons never oc- 
curs' M. 
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where r is the proper time of the particle. By sphtting the integral over x{t) on the right 
side of Eqn(1.8) into the product of a series of infinitesimal amplitudes corresponding 
to small segments, At, Gaussian integration may be performed over the intermediate 
space-time points. Finally, integrating over the proper time r, the analytical form of the 
propagator is found to be [0: 

K{xf, Xi] m) ^ ( ) Hi'\ims) (1.9) 

whereQ: 

S = ^{Xf - XiY 

and iff ^ is a first order Hankel function of the second kind, in agreement with Ref. H. 

In the asymptotic region where s ^ m~^, or for the propagation of on-shell parti- 
cles 0, the Hankel function reduces to an exponential and yields the configuration space 
propagator ~ exp(— zms) of Eqn(2.11) below. It is also shown in Ref. that energy and 
momentum is exactly conserved in the interactions and decays of all such 'asymptotically 
propagating' particles. The use of quasi-classical particle trajectories and the require- 
ment of exact energy-momentum conservation are crucial ingredients of the calculations 
presented below. 

The structure of the paper is as follows. In the following Section the case of neutrino 
or muon oscillations following pion decay at rest is treated. Full account is taken of the 
momentum wave-packets of the progagating neutrinos and muons resulting from the Breit- 
Wigner amplitudes describing the distributions of the physical masses of the decaying 
pion and daughter muon. The corresponding oscillation damping corrections and phase 
shifts are found to be very small, indicating that the quasi-classical (constant velocity) 
approximation used to describe the neutrino and muon trajectories is a very good one. 
The incoherent effects, of random thermal motion of the source pion, and of finite source 
and detector sizes, on the oscillation probabilities and the oscillation phases, are also 
calculated. These corrections are found to be small in typical experiments, but much 
larger than those generated by the coherent momentum wave packets. In Section 3, 
formulae are derived to describe neutrino oscillations following muon decay at rest or 
the /3-decay of radioacive nuclei. These are written down by direct analogy with those 
derived in the previous Section for pion decay at rest. In Section 4, the case of neutrino 
and muon oscillations following pion decay in flight is treated. In this case the two- 
dimensional spatial geometry of the particle trajectories must be related to the decay 
kinematics of the production process. Due to the non-applicability of the ultrarelativistic 
approximation to the kinematics of the muon in the pion rest frame, the calculation, 
although straightforward, is rather tedious and lengthy for the case of muon oscillations, 
so the details are relegated to an appendix. Finally, in Section 5, the positive aspects 
and shortcomings of previous treatments in the literature of the quantum mechanics of 
neutrino and muon oscillations are discussed in comparison with the method and results 
of the present paper. 



The metric for four-vector products is time-like. 
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2 Neutrino and muon oscillations following pion de- 
cay at rest 



To understand clearly the different physical hypotheses and approximations underly- 
ing the calculation of the particle oscillation effects it is convenient to analyse a precise 
experiment. This ideal experiment is, however, very similar to LNSD 

MEN m 



17 



and KAR- 

except that neutrinos are produced from pion, rather than muon, decay at 



rest. 



The different space-time events that must be considered in order to construct the 
probability amplitudes for the case of neutrino oscillations following pion decay at rest 
are shown in Fig.l. A 7r+ passes through the counter Ca, where the time to is recorded, 
and comes to rest in a thin stopping target T (Fig. la)). For simphcity, the case of only two 
neutrino mass eigenstates z/i and z/2 of masses mi and m2 {mi > 1712) is considered. The 
pion at rest constitutes the initial state of the quantum mechanical probability amplitudes. 
The final state is an e~p system produced, at time Id-, via the process —>■ e~p at a 
distance L from the decaying 7r+ (Fig. Id)). Two different physical processes may produce 
the observed e~p final state, as shown in Fig. lb) and Ic), where the pion decays either at 
time ti into z/i or at time t2 into z/2. The probability amplitudes for these processes are, 
up to an arbitary normalisation constant: 

A = /<e-p|T|n.,><.>.>D(x,-x.,t^-t.,m.)W(l^,,),m,,r,)<..|.,^ 

< z/^/i+|T|7r+ > e-^(*^-*°)Z}(0, U - to, m^)BW{W^, m^, T^)dW^ii) z = 1,2 (2.1) 

Note that following the conventional '/z' (final, initial) ordering of the indices of matrix 
elements in quantum mechanics, the path amplitude is written from right to left in order 
of increasing time. This ensures also correct matching of 'bra' and 'kef symbols in the 
amplitudes. In Eqn(2.1), < e~p|T|ni/e >, < z/^/i'^|T|7r+ > are the invariant amplitudes 
of the z/e charged current scattering and pion decay processes , BW(W^(i),m^, F^) and 
BW(WT,,mT,, Ftt) are relativistic Breit-Wigner amplitudes, < UeWi > and < z/i|z/^ > are 
amplitudes describing the mixture of the flavour and mass neutrino eigenstates, and D is 
the Lorentz-invariant configuration space propagator |]^, |^ of the pion or neutrino. The 
pole masses and total decay widths of the pion and muon are denoted by m^, F^ and m^, 
F^ respectively. For simplicity, phase space factors accounting for different observed final 
states are omitted in Eqn(2.1) and subsequent formulae. 

Because the amplitudes and propagators in Eqn(2.1) are calculated using relativistic 
quantum field theory, and the neutrinos propagate over macroscopic distances, it is a 
good approximation, as already discussed in the previous section, to assume exact energy- 
momentum conservation in the pion decay process, and that the neutrinos are on their 
mass shells, i.e. pi = mf, where Pi is the neutrino energy-momentum four- vector. In 
these circumstances the neutrino propagators correspond to classical, rectilinear, particle 
trajectories. 

The pion and muon are unstable particles whose physical masses W^^ and VF^(j) differ 
from the pole masses and appearing in the Breit-Wigner amplitudes and covariant 
space-time propagators in Eqn(2.1). The neutrino momentum Pj will depend on these 
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physical masses according to the relation: 



2W^ 



[rrii 



(2.2) 



Note that, because the initial state pion is the same in the two path amplitudes in 
Eqn(2.1) does not depend on the neutrino mass index i. However, since the pion decays 
resulting in the production of ui and 1/2 are independent physical processes, the physical 
masses of the unobserved muons, i — 1,2, recoiling against the two neutrino 

mass eigenstates are not, in general, the same. In the following kinematical calculations 
sufficient accuracy is achieved by retaining only quadratic terms in the neutrino masses. 



and terms linear in the small quantities : 5^^ = PV^ — 
simplification of the relativistic Breit-Wigner amplitudes: 

Fm 

BW{W,m,r) = — — — - 

— m^ + imV 

Tm 

~ 6{2m + 6) + imV 

= BW{d,r) + 0{S^) 
Developing Eqn(2.2) up to first order in mf , 6i and 6-„ yields the relation : 



W^M(j) ~ "^M- "^^^^ allows 



(2.3) 



Pi = ^0 



1 



[mt 



m 



2^2 



+ 



5^ 



(m: 



m. 



25m„ 



m„ {ml 
+ m 



ml) 



where 



mt 



m^[m; 



mlf 



ml -ml 
2m7r 



29.8MeV 



(2.4) 



(2.5) 



The term ~ S^^rn^ which is also included in Eqn(2.4) gives a negligible 0{mf) contribution 
to the neutrino oscillation formula. For muon oscillations, however, it gives a term of 
O(m^) in the interference term, as discussed below. Similarly, the exact formula for the 
neutrino energy: 



Wt - W'^^^ + m 



in combination with Eqn(2.4) gives, for the neutrino velocity: 



El 



1 - 



2Pi 



1 - 



2(5. 



{ml + ml) 



+ 



4:6im„ 



mAmt 



mf 



mt. 



mf 



+ 0{mt,5l,5^) 



(2.6) 



(2.7) 



This formula will be used below to calculate the neutrino times-of-fiight: t 



1,2. 



The unitary transformation specifying the mixing of the fiavour and mass neutrino 
eigenstates is defined by a single real mixing angle 9: 



We > 



cos^li^i > +sin^|i/2 > 



luu > — —sm9\ui > + cos ^1^2 > 



(2.8) 
(2.9) 
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where the flavour eigenstates |a >, a ■ 
respect the orthonormahty conditions: 



Ve, and the mass eigenstates \a >, a = z/i, z/2 



< a\P >= 5q,/3, < a\b >= 6, 



ab 



(2.10) 



The parts of the amplitudes requiring the most careful discussion are the invariant 
space-time propagators D, as it is mainly their treatment that leads to the different 
result for the neutrino oscillation phase found in the present paper, as compared to those 
having previously appeared in the literature. In the limit of large time-like separations, 
the propagator may be written as [0, |[: 



D{Ax, At,m) 



m 



27TiJ{Aty - (Ax) 



exp[-imJ {Aty - {Axf 



(2.11) 



D is the amplitude for a particle, originally at a space-time point (xj, ti), to be found at 
{xf, tf) and Ax = Xj — Xi, At = tf — ti. In the following, according to the geometry of the 
experiment shown in Fig.l, only one spatial coordinate will be considered (Ax = Xf — Xj) 
and only the exponential factor in Eqn(2.11), containing the essential phase information 
for particle oscillations will be retained in the amplitudes. Solid angle correction factors, 
taken into account by the factor in large brackets in Eqn(2.11), are here neglected, but 
are easily included in the final oscillation formulae. Writing then 



D{Ax, At, m) ^ exp[— zmw (At)^ — (Ax)^] = exp[— zmAr] = exp[— zA^] 



(2.12) 



it can be seen that the increment in phase of the propagator, A0, when the particle 
undergoes the space-time displacement (Ax, At) is a Lorentz invariant quantity equal 
to the product of the particle mass and the increment, Ar, of proper time. Using the 
relativistic time dilatation formula: 



At = 7Ar 



E 



At 



m 



and also the relation, corresponding to a classical, rectilinear, particle trajectory: 

L E 



At 



V 



(2.13) 



(2.14) 



gives, for the phase increments corresponding to the paths of the neutrinos and the pion 
in Fig.l: 



AC 



mjAr, 



^At. = 
5^ (ml + ml) 25im^ 



+ 



m^{ti - to) = m^{tD - to) - 



(2.15) 



m^ito - to) - m^L \ 1 + 



ml 
2PI 



^ _ 25^ [mj + m^) ^ 4(5,m^ 



{ml — mi) m^ 



(2.16) 
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Making the substitution U — to to — to — L/vi in the exponential damping factor due 
to the pion hfetime in Eqn(2.1) and using Eqns(2.3),(2.12),(2.15) and (2.16) Eqns(2.1) 
may be written as: 



f - r 



2 (S^ + i^) 



mj / m 



Po V2Po 



^-1 L 



2 (Si + i'-f) 
dSi i = 1,2 (2.17) 



where 



and 



bo = mn{L -to + to 



0£i 



(m^ — m^)^ 



2m^(m^ + mf^YL 



\m. 



ml + 



(2.18) 
(2.19) 

(2.20) 
(2.21) 



To perform the integral over Si in Eqn(2.17) it is convenient to approximate the mod- 
ulus squared of the Breit-Wigner amplitude by a Gaussian, via the substitution: 



r 

2 



r / 6 



5 + 



2 \S'' + 



r2 



-(<5-z-)exp -— 



(2.22) 



where the width of the Gaussian is chosen so that it has approximately the same full 
width at half maximum, F, as the Breit-Wigner function. After the substitution (2.22), 
the integral over 5i in Eqn(2.17) is easily evaluated by a change of variable to 'complete 
the square' in the argument of the exponential, with the result: 



9 /"oo 

= f/ 



i^)) exp l^-^t + ^oii^i^ d5i = ^yi^r/. exp 



12 



(2.23) 

Since Tj, = 2.53 x 10~^'' MeV, the term dependent on in Eqn(2.19) may be neglected, 
so that ctj may be taken as a real number with the numerical value: 



ai 



3.1 X 10''' ( L{m) MeV-' 



(2.24) 



For typical physically interesting values (see below) of mj = 1 eV and L = 30 m, a, takes 
the value 0.48 MeV-^ so that 



aiV^ = 0.48 X 3.00 x 10"^^ = 1.44 x 10"^^ 



Then, to very good accuracy, Ii — I2 — — iy 7r/3F^, independently of the neutrino mass. 
It follows that for neutrino oscillations, the muon mass dependence of the amplitudes may 
be neglected for any physically interesting values of rrii and L. 
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Prom Eqns (2.17) and (2.23) the probability to observe the reaction — > e p at 
distance L from the pion decay point and at time to is: 

P{e-p\L,tD) = \Ai + A2\^ 



3 



^1 < e~p\T\nu^ > 1^1 < i/^/i+|T|7r+ > 



2 



xsin^^cos^^e-^^^*^-*")^ (2.25) 



(tf'+t{') 
-2e " 2 _Re exp t 



^-l)^ + K(l)-«42)]5. 



The time dependent exponential factors in the curly brackets of Eqn(2.25) are easily 
understood. If mi > m2 then ti > ^2 ■ This implies that the neutrino of mass mi 
results from an earlier decay than the neutrino of mass m2, in order to be detected at the 
same time. Because of the exponential decrease with time of the pion decay amplitude, 
the contribution to the probability of the squared amplitude for the neutrino of mass 
mi is larger. The interference term resulting from the product of the decay amplitudes 
of the two neutrinos of different mass, has an exponential factor that is the harmonic 
mean of those of the squared amplitudes for each neutrino mass eigenstate, and so is also 
suppressed relative to the squared amplitude for the neutrino of mass mi. The integral 
over the physical pion mass is readily performed by replacing the Breit-Wigner function 



by a Gaussian as in Eqn(2.22). This leads to an overall multiplicative constant ^ti/?,Vt^ 
and a factor: 

exp[-(a,(l) - a,{2)fTl/l2] 

multiplying the interference term. For Am^ = (leV)^ and L — 30m the numerical value 
of this factor is exp(— 1.3 x 10~^^). This tiny correction is neglected in the following 
equations. 

Integrating over to gives the average probability to observe the e~p final state at 
distance L: 

P{e-p\L) = < e-p\T\nVe > Tl < i^^/x+|r|7r+ > l^sin^^cos^^ 

rVrtjAm^ 2m^m^Am^ | 

Where all kinematical quantities are expressed in terms of Am^, m^ and m^. Note that 
the minimum value of t^i is to + t\i to + 12 and to + ti for the squared amplitude terms of 
neutrinos of mass mi, m2 and the interference term, respectively. On integrating over to, 
the squared amplitude terms give equal contributions, the larger amplitude for mass mi 
being exactly compensated by a smaller range of integration. The exponential damping 
factor in the interference term in Eqn(2.26) is derived using the relations: 
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and 



2 

mf 



to obtain 



''^(''^ = '^-^ + 0^'^') ^ = 1,2 (2.28) 

fi fi _ {ml - ml)L _ Am^L 

The damping factor arises because the difference in the times-of-fiight of the two neutrino 
paths is hmited by the pion hfetime. It will be seen below, however, that for distances L 
of practical interest for the observation of neutrino oscillations, the damping effect is tiny. 

The part of the oscillation phase in Eqn(2.25) originating from the neutrino propa- 
gators (the term associated with the '1' within the large curved brackets) differs by a 
factor two from the corresponding expression in the standard formula. The contribution 
to the oscillation phase of the propagator of the decaying pion ( the term associated with 
m7r/(2Po) in the large curved brackets of Eqn(2.24)) has not been taken into account in 
any published calculation known to the author of the present paper. The oscillation phase 
in Eqn(2.26) is 2m^/(m^ — m^) = 2.685 times larger than that given by the standard 
formula (1.1). For L = 30m, as in the LNSD experiment, the first oscillation maximum 
occurs for Am^ = 0A6{eVy. Denoting by (f)i2 the phase of the cosine interference term 
in Eqn(2.26), the pion lifetime damping factor can be written as: 

F'^ir^) = exp |^-r,^0rr) = exp(-1.58 x 10-1^^2^) (2-30) 

so the damping effect is vanishingly small when (j)'^2 — 1. 

The oscillation formula (2.26) is calculated on the assumption that the decaying pion 
is at rest at the precisely defined position Xi. In fact, the positive pion does not bind with 
the atoms of the target, but will rather undergo random thermal motion. This has three 
effects: an uncertainy in the value of Xi, a Doppler shift of the neutrino energy and a time 
dilatation correction correction factor of 1/7^ in the equation (2.16) for the pion phase 
increment. Assuming that the target is at room temperature (T= 270° K), the mean 
kinetic energy of 2>kT /2 correponds to a mean pion momentum of 2.6 x 10^^ MeV and 
a mean velocity of ~ 5.6 km/sec. The pion will move, in one mean lifetime (2.6 x 10~^ 
sec), a distance of 146//m. This is neghgible as compared to L (typically > 30m) and so 
Eqn(2.26) requires no modification to account for this effect. 

The correction factor due to the Doppler effect and time dilatation is readily calculated 
on the assumption of a Maxwell-Boltzmann distribution of the pion momentum: 



dN 2 ( Vl 
— ~p^exp 



c.p^exp(-5 (2.31) 



Here — yfl^Tm^ — 2.64 x 10 ^ MeV. Details of the calculation are given in Appendix 
A. The interference term in Eqn(2.26) is modified by a damping factor: 



F''{Dop) = exp <^ - 



/ p^Am^ 



\ 2m^Po 



m 








(2.32) 
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while the argument of the cosine term acquires an additional phase factor: 



(j)''{Dop) 



3 / \ ^ Am^ r3m. 



4 Vm^y Po \-2Po 



1 L 



(2.33) 



For = 1, F-'iDop) = 1 - 6.7 X 10"^° and (p^iDop) = 1.2 x 10"^ 

If the target in which the pion stops is of thickness £t, then the effect of different 
stopping points of the vr (assumed uniformly distributed) is to multiply the interference 
term in (2.26) by the factor: 



If the position of the neutrino interaction point within the target has an uncertainy of 
±££)/2 a similar correction factor is found, with the replacement ix in Eqn(2.34). 

The calculation of this correction factor is also described in Appendix A. 

If the target T in which the pion comes to rest (Fig. la)) is chosen to be sufficiently thin, 
the pion decay process may be detected by observing the recoil muon in the counter Ca 
at times ti or t2 (Fig. lb) or Ic)). A sufficiently accurate measurement of the times ti, t2 
and to would, in principle, enable separation of the different processes in Figs, lb) and Ic) 
by the observation of separated peaks in the time-of-fiight distribution aXt{ = t/? — ti and 
^2 = tD — t2. In this case the interference term in Eqn(2.26) vanishes as the two alternative 
space-time paths leading to the neutrino interaction are distinguishable. However, for 
L = 30m and Am^ = l(eV)^ the difference in the times of flight is only 5.6 x 10"^^ 
sec, more than ten orders of magnitude smaller than can be measured with existing 
techniques. As discussed in Section 5 below, the momentum smearing due to the Doppler 
effect at room temperature is some eleven orders of magnitude larger than the shift due 
to a neutrino mass difference with Am? = l(eV)^. Thus, even with infinitely good time 
resolution, separation of such neutrino mass eigenstates by time-of-flight is not possible. 

The ideal experiment, described above to study neutrino oscillations, is easily adapted 
to the case of oscillations in the decay probability of muons produced by charged pion 
decay at rest. As previously pointed out in Ref. 0], such oscillations will occur if neutrinos 
with different masses exist. As before, the pion stops in the target T at time to (s^e 
Fig. 2a)). At time ti the pion decays into ui and the corresponding recoil muon (/ii), 
whose passage is recorded in the counter Cb (Fig. 2b)). Similarly, a decay into z/2 and /i2 
may occur at time t2 (Fig.2c)). With a suitable choice of the times ti and t2, such that 
muons following the alternative paths both arrive at the same time t^j at the point Xf, 
interference occurs between the path amplitudes when muon decay occurs at the space- 
time point {xf, tf)) in the detector D (Fig. 2d)). The path amplitudes corresponding to 
muons recoiling against neutrinos of mass mi and m2 are: 



Ai = y < e+z/eZ/^|T|/i+ > e D{xf - Xi^to - ti,m^)BW{W^(i),m^,T^) < > 
< iy^fi^\T\n+ > e-^(*^-*")D(0, ti - to, m^)BW{W^, m^, T^)dW^(^^) z = 1,2 (2.35) 



The various factors in these equations are defined, mutadis mutandi, as in Eqn(2.1). 




(2.34) 
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With the same approximations, concerning the neutrino masses and the physical pion 
and muon masses, as those made above, the velocity of the muon recoiling against the 
neutrino mass eigenstate Vi is: 



^0 



1 - 



+ 



{ml — rn?^'^{m'^ + m^) m: 



m. 



{mi - mty{ml - ml) 



(2.36) 



where 



^0 



m: 



m. 



(2.37) 



ml + 

Comparing with Eqn(2.7), it can be seen that for the muon case, unlike that where 
neutrino interactions are observed, there are pion and muon mass dependent correction 
terms that are independent of the neutrino masses, implying a velocity smearing effect 
due to the physical pion and muon masses that is ~ m^/mf larger than for the case of 
neutrino oscillations. 



The phase increments corresponding to the paths of the muons and the pion 
are, using (2.4)0 and (2.12)-(2.14): 



in Fig. 2 



A01 



mlL 



Po 

2 



1 + 



{ml + ml) 



2m^Pi 



6^ml {ml + ml) 



m 



{ml - m2)2 



{ml — ml 



+ 



mi — m^. 



(2.38) 



1^^^ = m^{ti - to) = m^{tD - to) - 



where 



mT,{t 



4(5, 



D 



-m^m 



m^^L 



1 + 



Amfmlmf^ 



imt 



m'j^)'^{ml + m^) 



mi — m^ 



^ 4:5im^ml ^ S^^^m^m^m^ (3m2 - m^ 



mi — mi 



{mi 



m 



m^i 



m^ 



ml + 
2m-^ 



(2.39) 
(2.40) 



Using Eqns(2.12),(2.38) and (2.39) to re- write the space-time propagators in Eqn(2.35), 
as well as Eqn(2.3) for the Breit-Wigner amplitudes gives: 



< e z/eZ/^|T|;U > e < z/j|z/^ >< i/^j^^ \T\7t > 



{Sn + 



exp t 



^Note that, in the pion rest frame Pi — 




2 {6, + i\) 2 

d5i i = l,2 (2.41) 
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where 



= 



^0 



-t 



D 



mt 



m. 



2 ml 



m. 



<(0 



2mlL 



ml -ml 



1 + 



m^(5m^ — lOm^m 



mlmf^ 



ml) 



and 



{mi - ml) {ml - ml)^ 
, / 1 



ml - ml 



t 



Hii) 



2m^Pl 



(2.42) 
(2.43) 
(2.44) 

(2.45) 



Making the substitution (2.22) and performing the integral over Si according to Eqn(2.23), 
the following formula is found for the probability for muon decay at distance L and time 
to. 

\A1 + A^\^ 

A'l2 1 ^ + — iT-il 4 



P{e'^UeUij_\L,tD) 



> 



+ \ |2 — r,r(t£)— to)_ 



m + ? ) 



sin^ ^e^"*''(i) + cos^ 0e 



m(2) 



-2sin6'cos6'e'^2"(*i(i)+*^(2)) 



Re exp i 



m^Am"^ 


2Pi 



L+K(l)-<(2)]5. 



(2.46) 



The muon path difference yields the term associated with Eq / mT^ in the interference phase 
in Eqn(2.46) while the pion path is associated with '1' in the large round brackets. The 
numerical value of the damping factor: 



exp[- 



{a^T, 



Ml2 



6 " 3 ' 

resulting from the integral over the physical muon mass is, for L = 30m, 0.774, so, unlike 
for the case of neutrino oscillations, the correction is by no means negligible. This is 
because, in the muon oscillation case, the leading term of is not proportional to the 
neutrino mass squared. The non-leading terms proportional to mf have been neglected 
in Eqn(2.43). This correction however effects only the overall normalisation of the oscil- 
lation formula, not the functional dependence on L arising from the interference term. 
Integrating over 6^, using Eqns(2.22) and (2.23), as well as over to, gives the probability 
of muon decay at distance L from the production point, where all kinematical quantities 
are expressed in terms of Am^, and m^: 



P{e+v,V^\L) 



ij. 

3a/3 



exp 



mt 



< e'^i'e^n\T\iJ' > I exp 



3 m^ — ml 



2T^m^m^{ml - ml) 



{ml + m2)3 



. 2T^mlmlAm^ 2m^mlAm^ 
sm 20 exp[ — L\ cos ■ 



{ml — ml)^ 



[mt 



m2)2 



< u^pi+\T\n+ > 



(2.47) 
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In this expression the correction due to the damping factor of the interference term: 



exp[-«(l) - <(2))2r^/12] 



arising from the integral over the physical pion mass has been neglected. For Am^ = 

(leV)^ and L — 30m the numerical value of this factor is exp(— 3.9 x 10^'^^). Denoting by 
(f)i2 the argument of the cosine in Eqn(2.47), the exponential damping factor due to the 
pion lifetime may be written as: 



exp 



{ml 



m 



I) 



(2.48) 



— 1, Fi^(Tt,) — 1 = 4.4 X 10 so, as in the neutrino oscillation case, the pion 
lifetime damping of the interference term is very small. 



For 0f 



Corrections due to time dilatation and the Doppler effect are calculated in a similar 
way to the neutrino oscillation case with the results (see Appendix A): 



F^{Dop) = exp 



. 2777,P3 



CO 










2 







and 



3 




V 2 777^A777^ 




2 




) Pi 


2r777r 



(2.49) 



(2.50) 



As for neutrino oscillations, the corresponding corrections are very small for oscillation 
phases of order unity. 

The phase of the cosine in the interference term is the same in neutrino and muon 
oscillations, as can be seen by comparing Eqns(2.26) and (2.47). It follows that the target 
or detector size correction (Eqn(2.34)) is the same in both cases. 

Neutrino and muon oscillations from pion decay at rest have an identical oscillation 
phase for given values of Ar77^ and L. In view of the much larger event rate that is possible, 
it is clearly very advantageous in this case to observe muons rather than neutrinos. In 
fact, it it not necessary to observe muon decay, as in the example discussed above. The 
oscillation formula applies equally well if the muons are observed at the distance L using 
any high efficiency detector such as, for example, a scintillation counter. In contrast, the 
rate of neutrino oscillation events is severely limited by the very small neutrino interaction 
cross section. 



3 Neutrino oscillations following muon decay at rest 
or beta-decays in nuclear reactors 

The formula describing neutrino oscillations 17^ — * z7g following the decay at rest of 
a /i"*", ^ e^Uf.Vfj_ is easily derived from the similar formula where a is produced 
by decay at rest (2.25). Because the neutrino momentum spectrum is continous, 
smearing effects due to the finite muon lifetime may be neglected from the outset. The 
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phase increment associated with the neutrino path is then given by Eqn(2.15) with the 
replacements Pq —>■ Py and 5T^,5i 0, where Py is the antineutrino momentum. The 
phase increment of the decaying muon is given by the same replacements in Eqn(2.16) 
and, in addition, and T^^ T^. The formula, analogous to Eqn(2.26), for the 

time-averaged probability to detect the process VeP —>■ e^n at a distance L from the muon 
decay point is then: 



P(e+n,/i|L) 



< e^n\T\pU(, > PI < Uf,i'^j.e'^\T\ij,'^ > 



,2 



X \ 1 -exp[ -^^L] cos 



Pu \2Ptj 



2 sin^ e cos^ e 



(3.1) 



The standard neutrino oscillation formula, hitherto used in the analysis of all exper- 
iments, has instead the expression Am'^L/{2Pij) for the argument of the cosine term in 
Eqn(3.1). Denoting my Am| the value of Am^ obtained using the standard formula, and 
Am'jpp that obtained using the Feynman Path (FP) formula (3.1) then: 

Am|p = (3.2) 
p- 

For a typical value of P-^ of 45 MeV, Eqn(3.2) implies that Arn^pp ~ 2.9Am|. Thus 



the oscillation signal from /i"*" decays at rest reported by the LNSD Collaboration 
corresponding to Am| ^ 0.5 (eV)^ for sin^2^^ ^ 0.02 implies Am^pp ~ 1.5 (eV)^ for a 
similar mixing angle. 

For the case of /?-decay: 

N{A, Z) N{A, Z + l)e-Ve 

771^^ in the first line of Eqn(2.16) is replaced by Ej^, the total energy release in the /5-decay 
process: 

Ep = MNiA,Z)-MNiA,Z+l) (3.3) 

where Mn{A, Z) and Mn{A, Z+1) are the masses of the parent and daughter nuclei. That 
the phase advance of an unstable state, over a time. At, is given by exp{—iE* At) where E* 
is the excitation energy of the state, is readily shown by the application of time-dependent 
perturbation theory to the Schrodinger equation [|T^. A more intuitive derivation of this 
result has been given in Ref. |2^. In the present case, E* = Ep. Omitting the lifetime 



damping correction, which is about eight orders of magnitude smaller than for pion decay, 
given a typical /3-decay hfetime of a few seconds, the time-averaged probabilty to detect 
the z/g via the process V^p — > e^n, at distance L from the decay point is given by the 
formula, derived in a similar way to Eqns(2.26) and (3.1): 



P{e+n,p\L) 



< e+n\T\pVe > PI < e-VeN{A, Z + l)|r|A^(A, Z) > 



X I sin^ 9 + cos^ 9 + 2 sin^ 9 cos^ 9 cos 



Ptj \2Py 



(3.4) 
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where F^^ = Tp the hfetime of the unstable nucleus N[A, Z). Until now, all experiments 
have used the standard expression t^n^Lji^Py) for the neutrino oscillation phase. The 
values of ISw? found should be scaled by the factor {Ep/ Py — 2)~^ ^ suitably averaged over 
Pu, to obtain the Am^ given by the Feynman Path formula (3.4). 



Neutrino and muon oscillations following pion de- 
cay in flight 



In this Section, the decays in flight of a vr"*" beam with energy E-,^ ^ into ^ are 
considered. As the analysis of the effects of the physical pion and muon masses have been 
shown above to give neghgible corrections to the L dependence of the oscillation formulae, 
for the case of decays at rest, such effects will be neglected in this discussion of in-flight 
decays. The overall structure of the path amplitudes for neutrinos and muons is the same 
as for decays at rest (see Eqns(2.1) and (2.35)). However, for in-fiight decays, in order to 
calculate the interfering paths originating at different and terminating at common space- 
time points, the two-dimensional spatial geometry of the problem must be properly taken 
into account. 

In Fig. 3 a pion decays at A into the 1 mass eigenstate, the neutrino being emitted 
at an angle 9i in the lab system relative to the pion flight direction. If mi > 1712 a later 
pion decay into the 2 mass eigenstate at the angle 9i + 5 may give a path such that both 
eigenstates arrive at the point B at the same time. A neutrino interaction u^n e^p, or 
i^fi^ — > e~p occuring at this space-time point will than be sensitive to interference between 
amplitudes corresponding to the paths AB and ACB. The geometry of the triangle ABC 
and the condition that the 1 and 2 neutrino mass eigenstates arrive at B at the same time 
gives the following condition on their velocities: 



Vii^u) sin6'2 t'i(z/) sin(^^2 — 6*1) 
V2{i') sin^^i sin^^i 



(4.1) 



Expanding to first order in the small quantity S = 62 — Oi, rearranging and neglecting 
terms of 0{m^) gives: 



V2{u) - vi{u) 



2E^ sin 9i 



1 — i'tt cos 9i 



where the relation: 



Vi{iy) = 1 



2E! 



has been used. Rearranging Eqn(4.2): 

Am^ 



2E?. 



sin 61 
1 — cos 9i 



(4.2) 



(4.3) 



(4.4) 



The difference in phase of the neutrino paths AB and CB is (see Eqn(2.15)): 

,^ mlAB mlCB 4 4, 



Pi 



(4.5) 
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Since the angle 5 is ~ Am^, the difference between AB and CD is of the same order, and 
so: 

'i'''^^-^rf^ + 0{m\,mt) (4.6) 

cos Uitji, 

where Pi P2 — Ey, the measured neutrino energy. From the geometry of the triangle 
ABC: 

AC AB L 



sin (5 sin^^2 cos^i sin(^i + (5) 
So, to first order in 5, and using Eqn(4.4): 



(4.7) 



j^Q ^ = — - f4 8) 

cos ^1 sin ^1 2£'2cos^i (1 — -UTrCOS^i) 

and 

Ax^ Am^L 1 
Ai^r = - = (4.9) 

Vtt 2£'^ COS^^i (1 — fjr COS^l) 

Eqns(4.8) and (4.9) give, for the phase increment of the pion path: 

Am^ EL (1 — v"^) 

Ar = m,(T2 - Ti) = m^Ar = E^At, - p^Ax, = \ — ^ ^ (4.10) 

2E^ cos &i (1 — cos Ui) 

In Eqn(4.10), the Lorentz invariant character of the propagator phase is used. Setting 
cos^i = 1 and = 0, the pion phase increment of Eqn(2.25) is recovered. Eqns(4.6) and 
(4.10) give, for the total phase difi^erence of the paths AB, ACB: 



Am^L 
cosOiEy 



1 - 



2Ey (1 -t;^ cos 6*1) 



(4.11) 



Using the expressions, valid in the ultra-relativistic (UR) limit where — 1 : 



771 2 1 

1-t;^cos^i = -^7 — (4.12) 

£;2(l + cose*) ^ ^ 

and 

E^{ml-ml) 

Eu= 2m=^ (1 + cosO (4.13) 

where 91 is the angle between the directions of the pion and neutrino momentum vectors 
in the pion rest frame, Eqn(4.11) may be rewritten as: 

Am^ ml 
COS uihii, [m^ — mj^) 

For ^1 = the oscillation phase is the same as for pion decay at rest (see Eqn(2.26)) 
since in the latter case, E^ c:^ Pq = (m^ — m^)/(2m^). Using Eqn(4.14), the probability 
to observe a neutrino interaction, at point B, produced by the decay product of a pion 
decay occuring within a region of length, lo^c, (-C L) centered at the point A, in a beam 
of energy Ej^, is given by a formula analagous to Eqn(2.26): 

P{e-p\L,0i) = !^££^^| < e-p\T\niye > H < i^^/^+|T|7r+ > l^sin^^cos^^ 

En 

{ nilAm'^L ] 
xU-cos^^^ — ^^ttt; 7r> 4.15 
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As in the case of pion decay at rest, Eqn(2.26), the oscillation phase differs by the factor 
2m^/ (m^ — mfj) — 2.685 from that given by the standard formula. 

The derivation of the formula describing muon oscillations following pion decays in 
flight is very similar to that just given for neutrino oscillations. The condition on the 
velocities so that the muons recoiling against the different neutrino mass eigenstates arrive 
at the point B (see Fig. 3) at the same time, is given by a formula analagous to (4.2): 



vi{li)[vi{li) --t;^cosgi] 

Av{fi) = V2{fi) - viin) = — 5 

sm Oi 



(4.16) 



The formula relating the neutrino masses to the muon velocities is, however, more difficult 
to derive than the corresponding relation for neutrinos, (4.3), as the decay muons are 
not ultra-relativistic in the pion rest frame. The details of this calculation are given in 
Appendix B. The result is: 



2m^E, 



[ml + ml)E^^ 



Using Eqn(4.16) and the relation, vahd to ffist order in 5,: 

At = 



cos 9i sin 9i 



(4.17) 



(4.18) 



where At is the flight time of the pion from A to C in Fig. 3 (and also the difference in 
the times-of- flight of the muons recoiling against the two neutrino eigenstates) , the angle 
S may be eliminated to yield: 



At 



Mf')L 



Vl(/i) COS^l[vi(/i) — COS^i] 

Using now the kinematical relation (see Appendix B): 



Vi{^) — VttCOS^i = _ „ „ ^ I 1 



{ml + ml)E^^ 



and the expression for the phase difference of the paths AB and ACB: 

A0^^ - A0^^^ 



At(^-^' 



E. 



Ett 



(4.19) 



(4.20) 



(4.21) 



together with Eqn(4.19), it is found, taking the UR limit, where i'i(//),i',r — 1, that 

{mlE^ - mlE,)L 



'12 



2mf,Am'^ 



El{ml - mlf 



cos 9\ 



(4.22) 



The probability of detecting a muon decay at B is then: 



X exp 



En COS 6*1 



2vr? Avr? 

sin 2Q cos 



[mlE, - mlE,)L 
cos^i 



■(4.23) 
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where Idbc is defined in the same way as in Eqn(4.13). 



5 Discussion 



The quantum mechanics of neutrino oscillations has been surveyed in two recent review 
articles , where further extensive lists of references may be found. 



The standard derivation of the neutrino oscillation phase will first be considered, 
following the treatment of Ref. [0, but using the notation of the present paper. The 
calculation is performed in terms of the 'flavour eigenstates of the neutrino', \a >, which 
in the case of pion decay is, at t = 0, essentially the pure state a = Ufj, corresponding 
to the superposition of the mass eigenstates |z/i >, |z/2 > given in Eqn(2.9) above. This 
flavour eigenstate is assumed to evolve with laboratoy time, t, according to fixed energy 
solutions of the non-relativistic Schrodinger equation: 

\a,t>= -sin^e-*^i*|z/i,j9 > + cos^e'^-^^Vz, P > (5.1) 

where |z/j,p > are mass eigenstates of fixed momentum p, and Ei, E2 are the laboratory 
energies of the neutrino mass eigenstates. The amplitude for transition into the state 
\i^e-,P > at time t is then, using Eqns(2.8), (2.9): 

< z/e, p\a, t >= sin 6 cos 6 (-6"*^^* + 6"*^^*) (5.2) 

Assuming now that the neutrinos have the same momentum but different energies : 

2 

/ 777 

= vy+^ = P+ir + Oimf) (5.3) 
and using (5.2) and (5.3), the probability of the flavour state at time t is found to be: 



P(z/e,t) = I < z/e,p|a,t > I = 2 COS 61 sin'' 6' 1 + cos 



[m^ — mr, 



2p 



(5.4) 



Finally, since the velocity difference of the neutrino mass eigenstates is O(Am^), then, to 
the same order in the oscillation phase, the replacement t — > L can be made in Eqn(5.4) 
to yield the standard oscillation phase of Eqn(l.l). 

The following comments may be made on this derivation: 

(i) The time evolution of the neutrino mass eigenstates in Eqn(5.1) according to the 
Schrodinger equation yields a non-Lorentz-invariant phase ~ Et, to be compared 
with the Lorentz-invariant phase ~ m?t/E given in Eqn(2.12) above. Although the 
two expressions agree in the non-relativistic limit ~ m it is clearly inappropriate 
to use this limit for the description of neutrino oscillation experiments. It may be 
noted that the Lorentz-invariant phase is robust relative to different kinematical 
approximations. The same result is obtained to order for the phase of spatial 
oscillations independent of whether the neutrinos are assumed to have equal mo- 
menta or energies. This is not true in the non-relativistic limit. Assuming equal 
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momenta gives the standard result of Eqn(l.l), whereas the equal energy hypothe- 
sis results in a vanishing oscillation phase. A contrast may be noted here with the 
standard treatment of neutral kaon oscillations, which follows closely the derivation 
in Eqns(5.1) to (5.4) above, except that the particle phases are assumed to evolve 
with time according to exp[— zmr] where m is the particle mass and r is its proper 
time, in agreement with Eqn(2.12). 

(ii) As pointed out in Ref. the different neutrino mass eigenstates do not have 
equal momenta as assumed in Eqns(5.1) and (5.3). The approximation of assuming 
equal momenta might be justified if the fractional change in the momentum of 
the neutrino due to a non-vanishing mass were much less than that of the energy. 
However, in the case of pion decay as is readily shown from Eqns(2.4) and (2.6) 
above, the ratio of the fractional shift in momentum to that in energy is actually 
(m^ + m^)/(m^ — m^J = 3.67; so, in fact, the opposite is the case. 

(iii) To derive (5.2) from (5.1), it is necessary to introduce the concept of a 'flavour 
momentum eigenstate' which cannot be defined in a theoretically consistent man- 
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(iv) What are the physical meanings of t, p in Eqn(5.1)? In Eqn(5.2), it is assumed that 
the neutrino mass eigenstates are both produced, and both detected, at the same 
times. Thus both have the same time-of-fiight t. The momentum p cannot be the 
same for both eigenstates, as assumed in Eqn(5.3), if both energy and momentum 
are conserved in the decay process. For any given value of the laboratory time t 
the different neutrino mass eigenstates must be at different space-time positions 
because they have different velocities[^ if it is assumed that both mass eigenstates 
are produced at the same time. It then follows that the Ue part of the different mass 
eigenstates cannot be probed at, some space-time point, by a neutrino interaction, 
whereas the latter must clearly occur at some unique and specific point in every 
event. 

(v) The historical development of the calculation of the neutrino oscillation phase is 
of some interest. The first published prediction p4| actually obtained a phase a 
factor two larger than Eqn(l.l) i.e. in agreement with the covariant calculation of 
Ref. (before the introduction of 'wave packets') and with the contribution from 
neutrino propagation found in the present paper. This prediction was later used. 



for example, in Ref. |25|. The derivation sketched above, leading to the standard 



result of Eqn(l.l) was later given in Ref. p6| . A subsequent paper 

in by the 



authors of Ref. p5| , published shortly afterwards, cited both Ref. [p^ and Ref.p6 
but used now the prediction of the latter paper. No comment was made on the 
factor of two difference in the two calculations. In a later review article PB[ by 



the authors of Ref. p5|| a calculation similiar to that of Ref. was presented in 



detail. Subsequently, all neutrino oscillation experiments have been analysed on the 
assumption of the standard oscillation phase of Eqn(l.l). 



It may be thought that the kinematical and geometrical inconsistencies mentioned 
in points (ii) and (iv) above result from a too classical approach to the problem. After 

^"This is true not only in the case of energy-momentum conservation, but also if it is assumed, as in the 
derivation of the standard formula, that the neutrinos have the same momentum but different energies. 
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all, what does it mean, in quantum mechanics, to talk about the 'position' and 'veloc- 
ity' of a particle, in view of the Heisenberg Uncertainty Relations ||2^? This point will 



become clear later in the present discussion, but first, following the original suggestion 



of Ref. and, as done in almost all subsequent work on the quantum mechanics of 
neutrino oscillations, the 'wave packet' description of the neutrino mass eigenstates will 
be considered. In this approach, both the 'source' and also possibly the 'detector' in 
the neutrino oscillation experiment are described by coherent spatial wave packets. Here 
the 'source' wavepacket treatment in the covariant approach of Ref. [0] will be briefly 
sketched. After discussing the results obtained, and comparing them with those of the 
present paper, the general consistency of the wave packet approach with the fundamental 
quantum mechanical formula (1.2) will be examined. 

The basic idea of the wave packet approach of Ref. [0] is to replace the neutrino prop- 
agator (2.11) in the path amplitude by a four- dimensional convolution of the propagator 
with a 'source wavefunction' which, presumably, describes the space-time position of the 
decaying pion. For mathematical convenience this wavefunction is taken to have a Gaus- 
sian form with spatial and temporal widths and at respectively. Thus, the neutrino 
propagator D is replaced by D where: 

(5.5) 



D{xf - Xi,mj) = I d xD{xf - x,,mj)ilJin{x - XiJ) 
where xt and Xf are 4- vectors and 

^in{x - Xi,j) = Noexp 



-iPj ■ [X - Xi) - 



4a| 



(5.6) 



The integral in (5.5) was performed by the stationary phase method, yielding the result 
(up to multiplicative and particle flux factors), and here assuming, for simplicity, one 
dimensional spatial geometry: 



D{Ax, At,mj) = exp 



-i y-^'^t - Pj{Ax - VjAt) 



{Ax - VjAtf 
4(^FmP) 



(5.7) 



For the case of — > z/g oscillations, following 7r+ decay at rest, the probability to observe 
flavour z/g at time t and distance x is given by: 



P(z/e, As, At) = cos^ e sin^ -D{Ax, At, nii) + D{Ax, At, ms) 



(5. 



Performing the integral over At and making the ultra-relativistic approximation vi,V2 — 1 
yields finally, with Ax = L: 



P{ue,L) = 2cos2^sin2e(l + 
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(5.9) 



It can be seen that the oscillation phase in Eqn(5.9) is the same as standard one of 
Eqn(l.l). The exponential damping factors in Eqn(5.9) are the same as those originally 
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found in Ref. 0] for spatial wave packets {at = 0). Considering now only spatial wave 
packets and using the property axCp = 1/2 derived from the Fourier transform of a 
Gaussian, the two terms in the exponential damping factor may be written as: 



exp 



and 



exp 



(5.10) 



(5.11) 



The spatial damping factor is usually interpreted in terms of a 'coherence length" 



If 



L > 



2(1711 — mlYa^ 



Am^(m^ + m^) ^ 

then <^ 1 and the neutrino oscillation term is strongly suppressed. Eqn(5.12) expresses 
the condition that oscillations are only observed provided that the wave packets overlap. 
Since Av ~ Am^ the separation of the wave packets is ~ AvL ~ Am^L, so that Eqn(5.12) 
is equivalent to AvL ^ cr^ (no wave packet overlap). 



The damping factor Fp is typically interpreted |33| in terms of the 'Heisenberg Un- 
certainty Principle'. This factor is small, unless the difference in mass of the eigenstates 
is much less than 24 ^m^cip, so it is argued that only for wide momentum wave pack- 
ets can neutrino oscillations be observed, whereas in the contrary case, when the mass 
eigenstates are distinguishable, the interference effect vanishes. In the case of pion decay 
the difference in momentum of the two interfering mass eigenstates comes only from the 
6i term in Eqn(2.4), as 6.,^, being a property of the common initial state, is the same for 
both eigenstates. The neutrino momentum smearing in pion decay is then estimated from 
Eqn(2.4) as: 

2Pn5ima r„m,, m 
2m^ 



Pi - Pu 



mf 



1.14 X IQ-'W 



For Am? = l(eV)^ the value of Fp is found to be exp(— 500) ~ 7.1 x 10"^^^ giving complete 
suppression of neutrino oscillations. This prediction is in clear contradiction with the tiny 
damping corrections found in the path amplitude analysis in of Section 2 above. 

The physical relevance of the damping factors F^ and Fp will be reviewed after the 
examination of the physical meaning of the wave packet convolution formula (5.5) which 
now follows. 



The convolution formula (5.5) may be re-written so as to isolate the contribution from 
the region: 

, Ax'i , Ax\ 

x[-^<xi< x\ + 

where x'^ is an arbitary fixed value of the coordinate xi (say, the spatial position of source 
pion parallel to the direction of neutrino propagation): 



22 



D{xf - Xi,mi) = f Dipin / ' ^ dxi+ [ dxi d^x + Ax\ [ [Dipin]^^^^, d 

J J-oo Jx[-\ J 1 

= f Dij.^d^x + Ax['miJ~{x[) (5.13) 



Here the X2^ x^, x^ averaged quantity Dijjin{x[) is a function of {xf)k — {xi)k, = 2, 3, 4 
as well as — x[ and x[ — Similarly: 

D{xf - Xi,m2) = Di/jind^x + Ax'lDipinix'l) (5.14) 

J{-Ax'{) 

Now in calculating the oscillation probability (5.8), Eqns(5.13) and (5.14) will give a 
contribution to the interference term from the terms proportional to Ax[, Ax'{ of: 



/(Ax'i, Ax'i) = Ax[Ax'l2ReD^in{x[yD^in{xi) (5.15) 

This contribution corresponds to interference between a pion decaying at an arbitary 
point Xi = x'l with another decaying at another arbitary point This is in 

contradiction with Eqn(1.2), according to which, the initial quantum state must be the 
same for all interfering path amplitudes. Referring to Fig 1., it can be seen that the 
different path amplitudes correspond to two (classically alternative) decay histories of the 
same pion which must clearly be at some unique spatial position at time to- 'Each pion 
interferes only with itself. Interference between two different pions never occurs' ||16|| . 
Thus the convolution, at amplitude level, of the neutrino propagator with 'source' and/or 
'detector' wavefunctions does not respect Heisenberg and Feynman's fundamental law of 
quantum mechanics (1.2). 

The preceding discussion of the derivation of the standard formula for the oscillation 
phase (1.1) in terms of 'flavour eigenstates' revealed contradictions and inconsistencies if 
the neutrinos are assumed to follow classical space-time trajectories. The 'source wave 
packet' treatment gives the standard result for the oscillation phase and predicts that 
the interference term will be more or less damped depending on the widths in configura- 
tion and momentum space of the wave packets. However the convolution formula (5.5) 
is evidently at variance with Eqn(1.2) and the requirement that particles can interfere 
only 'with themselves'. So do wave packets actually play a role in the correct quantum 
mechanical description of neutrino oscillations, as suggested in Ref.[pU|? Are the packets 



actually constrained by the Heisenberg Uncertainty Relations? How do the properties 
of the wave packet effect the possiblity to observe neutrino oscillations? The answers to 
all these questions are contained in the results of the calculations presented in Section 2 
above. They are now reviewed, with special emphasis on the basic assumptions made and 
the physical interpretation of the equations. 

Referring again to Fig.l, in a) a single pion comes to rest in the stopping target T. 
The time of its passage is recorded by the counter Ca, which thus defines the initial state 
as a TT"*" at rest at time to. This pion, being an unstable particle, has a mass that 
is, in general, different from its most likely value which is the pole mass m^r. What are 
shown in Fig. lb) and Fig.lc) are two different classical histories of this very same pion. 
In b) it decays into the neutrino mass eigenstate ui at time ti, and in c) it decays into the 
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Source 


Am = leV 






Doppler Effect 




4.4x10-^6 


3.8xl0-^« 


3.3x10-1*^ 


1.9x10-5 



Table 1: Different contributions to neutrino momentum smearing in pion decay at rest 
(see text). 

neutrino mass eigenstatc z/2 at time t2- Because these arc independent classical histories, 
the physical masses VF^(2) of the muons recoiling against the mass eigenstates ui 

and 1^2 respectively, are, in general, not equal. Taking into account now exact energy- 
momentum conservation in the decay processes (appropriate because of the covariant 
formulation used throughout) the eigenstates z/i and z/2 will have momenta which depend 
on and and VF^ and VF^(2) respectively. These momenta are calculated in 

Eqn(2.4). The distributions of Wt^ and are determined by Breit-Wigner amplitudes, 
that are the Fourier transforms of the corresponding exponential decay laws, in terms of 
the decay widths and respectively. In accordance with Eqn(1.2), only the Breit- 
Wigner amplitudes corresponding to the physical muon masses are (coherently) integrated 
over at the amplitude level. The integral over W„ (a property of the initial state) is 
performed (incoherently) at the level of the transition probability. Because of the long 
lifetimes of the tt and the corresponding momentum wave packets are very narrow, 
so that all corrections resulting from integration over the resulting momentum spectra 
are found to give vanishingly small corrections. Indeed, as is evident from the discussion 
of Eqn(5.10) above, the width of the momentum wave packet is much smaller than the 
difference in the momenta of the eigenstates expected for experimentally interesting values 
of the neutrino mass difference (say, Am^ — l(eV)^). However, contrary to the prediction 
of Eqn(5.10), this does not at all prevent the observation of neutrino oscillations. This 
is because the oscillations result from interference between amplitudes corresponding to 
different propagation times, not momenta, of the neutrino mass eigenstates. Tabic 1 
shows contributions to the fractional smearing of the neutrino momentum from different 
sources: 

(1) neutrino mass difference Am^ = l(eV)^ (Eqn(2.4)) 

(2) coherent effect due to the physical muon mass (5^ = F^/2 in Eqn(2.4)) 

(3) incoherent effect due to the physical pion mass (5^ = F^/2 in Eqn(2.4)) 

(4) incoherent Doppler effect assuming — 2.6 x IQ-^MeV 

The pion mass effect is of the same order of magnitude as the neutrino mass shift. The 
muon mass effect is two orders of magnitude smaller, while the Doppler effect at room 
temperature gives a shift eleven orders of magnitude larger than a (leV)^ neutrino mass 
difference squared. 

According to the usual interpretation of the Heisenberg Uncertainty Principle, the 
neutrinos from pion decay, which, as has been shown above, correspond to very narrow 
momentum wave packets, would be expected to have a very large spatial uncertainty. 
Indeed, interpreting the width of the coherent momentum wave packet generated by the 
spread in according to the the momentum-space Uncertainty Relation ApAx — 1 
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gives Ax = 1.27 km. Does this accurately represent the knowledge of the position of a 
decay neutrino obtainable in the experiment shown in Fig.l? Without any experimental 
difficulty, the decay time of the pion can be measured with a precision of 10~^''sec, by 
detecting the decay muon. Thus, at any later time the distance[3 of the neutrino from 
the decay point is known with a precision of c x 10~^°cm = 3cm. This is a factor 4x10^ 
more precise than the 'uncertainty' given by the Heisenberg relation. It is clear that the 
experimental knowledge obtainable on the position of the neutrino is essentially classical, 
in agreement with the theoretical description in terms of classical particle trajectories 
in space-time. In this case the momentum-position Uncertainty Relation evidently does 
not reflect the possible experimental knowledge of the position and momentum of the 
neutrino. This is because it does not take into account the prior knowledge that the mass 
of the neutrino is much less than that of the pion or muon, so that its velocity is, with 
negligible uncertainty, c. In spite of this, a Heisenberg Uncertainty Relation is indeed 
respected in the pion decay process. The Breit-Wigner amplitude that determines the 
coherent spread of neutrino momentum is just the Fourier transform of the exponential 
decay law of the muon. The width parameter of the Breit-Wigner amplitude and the 
muon mean lifetime do indeed respect the energy-time Uncertainty Relation T^r^ = 1. 

It is then clear, from this careful analysis of neutrino oscillations following pion decay 
at rest, that, in contradiction to what has been almost universally assumed until now, 
the neutrinos are not described by a coherent spatial wave packet. There is a coherent 
momentum wave packet, but it is only a kinematical consequence of a Breit Wigner 
amplitude. If for mathematical convenience, the momentum wave packet is represented 
by a Gaussian, a conjugate (and spurious) Gaussian spatial wave packet will be generated 
by Fourier transformation. Indeed, in the majority of wave packet treatments that have 
appeared in the literature, Gaussian momentum and spatial wave packets related by a 
Fourier transform with widths satisfying the 'Uncertainty Relation' apax = 1/2 have 
been used. The above discussion shows that such a treatment does not correspond to the 
actual knowledge of the space-time events that constitute realistic neutrino oscillation 
experiments 

The limitation on the detection distance L, for observation of neutrino oscillations, 
given by the damping factor F'^ir^^) of Eqn(2.30) is easily understood in terms of the 
classical particle trajectories shown in Fig.l. For a given velocity difference, the time 
difference ^2 — ti becomes very large when both neutrinos, in the alternative classical 
histories, are required to arrive simultaneously at a far distant detector. Because of the 
flnite pion lifetime, however, the amplitude for pion decay at time t2 is smaller than that 
at ti by the factor exp[— (^2 — ti)r^/2]. Integrating over all decay times results in the 
exponential damping factor F^ITt^) of Eqn(2.30). It is clear that, contrary to the damping 
factor Fx of Eqn(5.11), the physical origin of the L dependent damping factor is quite 
unrelated to 'wavepacket overlap'. 

Here, 'distance' is defined without specifying tlie flight direction of the neutrino. Precise, simultane- 
ous, measurement of the flight direction of the recoil muon also determines, with a similar precision, the 
neutrino direction. Evidently the thickness of the stopping target may be chosen sufficiently small that 
its contribution to the neutrino position uncertainty is negligible. 

^^It follows from this that discussions of quantum mechanical coherence in neutrino oscillations based 
on the properties of Gaussian wavepackets, as in Ref. [|4|, do not address the actual physical basis of 
neutrino oscillation experiments. A similar remark applies to the more recent and very extensive study 
of Ref. ||. 
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For a given value of t2 — ti the coherent neutrino momentum spread originating in the 
Breit Wigner amphtude for PV^ produces a corresponding velocity smearing that reduces 
the number of possible classical trajectories arriving at the detection event. This effect is 
taken into account in the integral shown in Eqn(2.23). The effect is shown, in Section 2 
above, to be much smaller than the (already tiny) pion lifetime damping described above, 
and it is neglected in Eqn(2.26). 

Finally, as descibed after Eqn(2.25), the (incoherent) integration over the Breit Wigner 
function containing Wj^ gives an additional damping correction to the interference term 
that is also very tiny compared to that due to the pion lifetime. 

As promised above, a final remark is now made on the physical interpretation of 
the damping factors (5.10) and (5.11) that have often been derived and discussed in 
the literature. Although as pointed out above, and to be further discussed below, the 
derivation of these factors is incorrect, as it results from coherent integration over different 
spatial positions of the neutrino source, similar factors do appear when performing the 
integrations over the (unobserved) neutrino momentum distributions. is replaced by 
F'^ITtt) and Fp by the factor resulting from the coherent integration over the physical 
muon mass W^. The reason for the huge suppression factor predicted by Eqn(5.10), and 
the tiny one found in the path amplitude calculation, is that, in deriving (5.10) and (5.11), 
it is assumed that the neutrino eigenstates are both produced and detected at equal times. 
This will only be possible if both the hypothetical 'wave packet overlap' is appreciable 
(Eqn(5.11)) and also the momentum smearing is sufficiently large that the time-of-flight 
differences due to the different neutrino masses are washed out (Eqn(5.10). In the path 
amplitude calculation interference and hence oscillations are made possible by different 
decay times of the source pion, and the damping factors analagous to F^ and Fp turn out 
to give vanishingly small corrections to the oscillation term. 

It may be remarked that the physical interpretation of 'neutrino oscillations' provided 
by the path amplitude description is different from the conventional one in terms of 'flavour 
eigenstates'. In the latter the amplitudes of different flavours in the neutrino are supposed 
to vary harmonically as a function of time. In the amplitudes for the different physical 
processes in the path amplitude treatment there is, instead, no variation of lepton flavour 
in the propagating neutrinos. If the mass eigenstates are represented as superpositions of 
flavour eigenstates using the inverses of Eqns(2.8) and (2.9), there is evidently no temporal 
variation of the lepton flavour composition. Only in the detection process itself are the 
flavour eigenstates projected out, and the interference effect occurs that is described as 
'neutrino oscillations'. In the case of the observation of the recoil muons no such projection 
takes place, but exactly similar interference effects are predicted to occur. As previously 
emphasised [||, the 'flavour oscillations' of neutrinos, neutral kaons and b-mesons are just 
special examples of the universal phenomenon of quantum mechanical superposition that 
is the basis of Eqn(1.2). 

A more detailed critical review will now be made of previous treatments of the quantum 
mechanics of both neutrino and muon oscillations in the literature. 

By far the most widespread difference from the path amplitude treatment of the 
present paper is the non-respect of the basic quantum mechanical formula, (1.2), by the 
introduction of wave packets to describe 'source' and/or 'detector' particles[7,8, 10, 11,30- 
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46]. Since in any practical neutrino oscillation experiment a single initial or final quantum 
state, as specified by Eqn(1.2), is not defined, but rather sets of initial and final states 
/ = and F = Y^mfm determined by experimental conditions, Eqn(1.2) may be 

generalised to: 



FI 



EE 

m I 



ki k2 k„ 



to be contrasted with the formula used in the references cited above: 
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(5.16) 



(5.17) 



Here, ipi and ipf are 'source' and 'detector' wave packets respectively. In Eqn(5.17) the ini- 
tial states ii and final states fm of Eqn(5.16), that correspond to different spatial positions, 
and also, possibly, different kinematic properties, of the source particle or detection event, 
are convoluted with spatial and/or temporal 'wave packets', that, as discussed above, are 
devoid of physical significance. Use of Eqn(5.17) leads to the absurd prediction (see, for 
example, Eqn(5.15) above) that different initial or final state particles interfere with each 
other. A possible reason for the widespread use of Eqn(5.17) instead of (5.16) may be 
understood following a remark of the author of Ref. |^ concerning a 'paradox' of the 
complete quantum field theory calculation that takes into account, by a single invariant 
amphtude, production, propagation and detection of the neutrinos. It was noticed that, 
if the amplitude for the complete chain of processes is considered to correspond to one 
big Feynman diagram, then integration over the space-time coordinates of the initial and 
final states will reduce the exponential factors, containing the essential information on 
the interference phase, to energy-momentum conserving delta functions, and so no oscil- 
lations will be possible. A related remark was made by the authors of Ref. ||11| who stated 



that, as they were assuming exact energy-momentum conservation, the integration over 
the space-time coordinates could be omitted. They still, however, (quite inconsistently, in 
view of the previous remark) retained the exponential factors containing the interference 
phase information. These considerations indicate a general confusion between momen- 
tum space Feynman diagram calculations, where it is indeed legitimate to integrate, at 
the amplitude level, over the unobserved space-time positions of the initial and final state 
particles, and the case of neutrino oscillations, where the amplitude must be defined in 
configuration space. In the latter case, it is the unobserved momenta of the propagating 
particles that should be integrated over, as is done in Eqns(2.1) and (2.35) above, and 
not the space-time positions of the 'source' or 'detector' particles, as in Eqn(5.17). 

It is clear that exact energy-momentum conservation plays a crucial role in the path 
amplitude calculation. This is valid only in a fully covariant theory. Still, several authors, 
in spite of the ultra-relativistic nature of neutrinos, used a non-relativistic theory to 
describe the production, propagation and detection of neutrinos ^ As is 

well known, in such 'Old Fashioned Perturbation Theory' energy is not conserved 
at the level of propagators and so no precise analysis of the kinematics and the space- 
time configurations of the production and detection events, essential in the covariant path 
amplitude analysis, is possible. 

Even when, in some cases, the complete production, propagation and detection process 
of the neutrinos were described[35-37,39,40], equal neutrino energies [37,39,40], equal neu- 
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trino momenta [36] or either |^ were assumed, in contradiction with energy-momentum 



conservation and a consistent space-time description of the production and detection 
events. As follows directly from Eqn(5.3) (or the similar formula, for the neutrino mo- 
mentum, obtained by assuming equal neutrino energies), in all the above cited references, 
the standard formula (1.1) for the oscillation phase was obtained. 

An interesting discussion of the interplay between different kinematical assumptions 
(not respecting energy-momentum conservation) and the space-time description of the 
production and detection events was provided in Ref. [^. This treatment was based 
on the Lorentz-invariant propagator phase of Eqn(2.12). By assuming either equal mo- 
mentum or equal energy for the propagating neutrinos, but allowing different times of 
propagation for the two mass eigenstates, values of ^^'2 agreeing with Eqns(2.15) and 
(2.25) above were found, i.e. differing by a factor two from the standard formula. Al- 
teratively, assuming equal velocities, (and hence equal propagation times) the standard 
result (1.1) was obtained. In this latter case, however, the masses, momenta and energies 
of the neutrinos must be related, up to corrections of 0{mj), according to: 

Since the ratio of the neutrino masses may take, in general, any value, so must then the 
ratio of their momenta. For the case of neutrino production from pion decay at rest with 
mi, 7712 -C m^,m^ the relation (5.18) is clearly incompatible with Eqn(2.4) which gives: 

P2 {ml - m2)2 

On the other hand, Eqn(5.19) is clearly compatible (up to corrections of 0(Am^)) with 
the 'equal momentum' hypothesis. There is similar compatiblity with the 'equal energy' 
hypothesis. Even so, the authors of Ref. |^ recommended the use of the equal velocity 
hypothesis. With the hindsight provided by the path amplitude analysis, in which the 
two neutrino mass eigenstates do indeed have different propagation times, it can be seen 
that the kinematically consistent 'equal momentum' and 'equal energy' choices are good 
approximations and the neutrino oscillation phase, resulting from the propagation of the 
neutrinos alone, is indeed a factor two larger than the prediction of the standard formula. 



Only the author of Ref. included the propagator of the decaying pion in the 
complete production-propagation-detection amplitude; compare Eqn(8) of Ref. with 
Eqn(2.1) above. However, no detailed space-time analysis of production and detection 
events was performed. Square spatial wave packets for the 'source' and 'detector' were 
convoluted at amplitude level as in Eqn(5.17). As the neutrinos were assumed to have 
equal energies the standard result for the oscillation phase was obtained. 



Although, following Ref. |^ , most recent studies of the quantum mechanics of neu- 
trino oscillations have considered the complete production-propagation-detection process, 
some authors still use, in spite of the criticisms of Ref. P3[, the 'flavour eigenstate' de- 



scription 0, ^ In the last two of these references a 'quantum field theory' approach 
ia adopted, leading to an 'exact' oscillation formula that does not make use of the 
usual ultra-relativistic approximation. This work did not include either exact kinematics, 
or an analysis of the space-time structure of production and detection events. Also the 
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effect of tlie propagator of tlie source particle was not considered. In Ref. |^ tlie equal 
energy liypotliesis was used and in Ref. [|^, ^ tlie equal momentum hypothesis. As a 
consequence, in all three cases the standard result was found for the oscillation phase in 
the ultra-relativistic limit. 

Correlated production and detection of neutrinos and muons produced in pion decay 
were considered in Ref. [Q. The introduction to this paper contains a valuable discussion 
of the universality of the 'particle oscillation' phenomenon. It is pointed out that this is 
a consequence of the general principle of amplitude superposition in quantum mechanics, 
and so is not a special property of the K° — K ,3*^ — 6 and neutrino systems which are 
usually discussed in this context. This paper used a covariant formalism that employed the 
'energy representation' of the space-time propagator. In the introduction, the important 
difference between Eqns(5.16) and (5.17) was also touched upon: 

' The reader will agree that one should not integrate over space if one is interested in 
spatial interference (or oscillation).' 

Even so, in the amplitude for the correlated detection of the muon and neutrino 
(Eqn(2.10) of Ref. ^) not only are the space-time positions of the production points of 
the neutrino and muon integrated over, but they are assumed to be at different space- 
time points. The propagator of the decaying pion is not included in the amplitude, and 
although exact energy-momentum conservation is imposed, no space-time analysis of the 
production and decay points is performed. Correlated spatial oscillations of neutrinos 
and muons are predicted, though with interference phases different from the results of 
both the present paper and the standard formula. Integrating over the spatial position 
of either the detected neutrino or muon destroys the interference term, so it is predicted 
that only correlated neutrino-muon oscillations may be observed, following pion decay, 
in contradiction with the results of the present paper, and, (in the case of neutrino os- 
cillations) with all previous studies^. Pion and muon lifetime effects were mentioned in 
Ref. but neither the role of the pion lifetime in enabling different propagation times 
for the neutrinos nor the momentum smearing, induced by the Fourier-transform-related 
Breit Wigner amplitudes, were discussed. 

The claim of Ref. M) that correlated neutrino-muon oscillations should be observable 



in pion decay was questioned in Ref. The authors of the latter paper attempted to 
draw conclusions on the possibility, or otherwise, of particle oscillations by using 'plane 
waves', i.e. energy-momentum eigenfunctions. As is well known, such wavefunctions are 
not square integrable, and so can yield no spatial information. The probability to find a 
particle described by such a wave function in any finite spatial volume is zero. Due to the 
omission of the (infinite) normalisation constants of the wavefunctions many of the equa- 
tions in Ref. |Tl| are, as previously pointed out |]50|, dimensionally incorrect. Momentum 



wavepackets for the decaying pion convoluted at amplitude level as in Eqn(5.17) were 
also discussed in Ref. |Tl|. Although exact energy-momentum conservation constraints 
were used, it was assumed, as in Ref. [||, that the muons and the different neutrino mass 
eigenstates are both produced and detected at common points (Eqn(35) of Ref. [|n|]). The 
latter assumption implies equal velocities, yielding the standard neutrino oscillation phase 
as well as the inconsistent kinematical relation (5.18). The authors of Ref. |TI| concluded 



-'^^ Although this conclusion follows from Eqn(3.22) of Ref. M, it was not pointed out in the paper. 
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that: 



(a) correlated n — u oscillations of the type discussed in Ref. could be observed, 
though with different oscillation phases. 

(b) oscillations would not be observed if only the muon is detected 

(c) neutrino oscillations can be observed even if the muon is not detected. 



Conclusion (b) is a correct consequence of the (incorrect) assumption that the muons 
recoiling against the different neutrino mass eigenstates have the same velocity. As both 
muons have the same mass they will have equal proper time increments. So according to 
Eqn(2.12) the phase increments will also be equal and the interference term will vanish. 
The conclusion (c) is in contradiction to the prediction of Eqn(3.22) of Ref. p[, according 
to which, no neutrino oscillations can be observed when the decay position of muon is 
integrated over. The path amplitude calculation of the present paper shows that conclu- 
sion (b) is no longer valid when the different possible times of propagation of the recoiling 
muons are taken into account. 



Observation of neutrino oscillations following pion decay, using a covariant formal- 
ism( Schwinger's parametric integral representation of the space-time propagator) was 
considered in Ref. |^. Exact energy- momentum conservation was imposed, but integra- 



tion over the pion spatial position at amplitude level, as in Eqn(5.17), was done and no 
account was taken of the contribution of the pion space-time propagator to the oscilla- 
tion phase. As the different eigenstates were assumed to be be produced or detected at 
identical space-time points equal propagation velocities were implicitly assumed, so that 



just as in Refs[g, g, |ri[], where the same assumption was made, the standard neutrino 
oscillation phase was obtained. In the conclusion of this paper the almost classical nature 
of the space-time trajectories followed by the neutrinos was stressed, although this was 
not taken to its logical conclusion in the previous discussion, e.g. the kinematical incon- 
sistency of the equal velocity hypothesis that requires the evidently impossible condition 
(5.18) to be satisfied. 



In a recent paper , the standard neutrino oscillation formula with oscillation phase 



given by Eqn(l.l) was compared with a neutrino decoherence model. In order to take into 
account incertainties in the position of the source and the neutrino energy, an average was 
made over the quantity L/AE^, assuming that it is distributed according to a Gaussian 
with mean value i and width a. The average was performed in an incoherent manner. 
Thus the calculation is closely analagous to those for the effects of target or detector length 
or of thermal motion of the neutrino source, presented in the Appendix A of the present 
paper. Perhaps uniquely then, in the published literature, in Ref. |51| the effects of source 
position and motion are taken into account correctly, according to Eqn(5.16) instead of 
Eqn(5.17). However, the source of the neutrino energy uncertainty is not specified. In as 
far as it is generated from source motion the calculation is, in principle, correct. There 
is however also the (typically much smaller, see Table 1 above, for the case of pion decay 
at rest) coherent contribution originating from the variation in the physical masses of the 
unstable particles produced in association with the neutrino, as discussed in detail above. 
It was concluded in Ref. that the Gaussian averaging procedure used gave equivalent 



results to the decoherence model for a suitable choice of parameters. 
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It is clearly of great interest to apply the calculational method developed in the present 
paper to the case of neutral kaon and b-meson oscillations. Indeed the use of the invariant 



path amplitude formalism has previously been recommended for experiments involv- 
ing correlated pairs of neutral kaons. Here, just a few remarks will be made on the main 
differences to be expected from the case of neutrino or muon oscillations. A more detailed 
treatment will be presented elsewhere [p3 |. 



In the case of neutrino and muon oscillations, the interference effect is possible as 
the different neutrino eigenstates can be produced at different times. This is because the 
decay lifetimes of all interesting sources (pions, muons, /3-decaying nuclei) are much longer 
than the time difference beween the paths corresponding to the interfering amplitudes. To 
see if a similar situation holds in the case of Ks — Kl oscillations, three specific examples 
will be considered with widely differing momenta of the neutral kaons: 

(I) ^ KsKl 

(II) n-p AK° at = 2 GeV 
(III) n-p AKO at = 10 GeV. 



These correspond to neutral kaon momenta of 108 MeV, 750 MeV and 5 GeV respectively. 
In each case the time difference {Atx) of production of Kg and Kl mesons, in order that 
they arrive at the same time at a point distant cjkTs (where is the usual relativistic 
parameter) from the source is calculated. Exact relativistic kinematics is assumed and 
only leading terms in the mass difference ArriK = ms — rriL are retained. Taking the value 



of ArriK and the various particle masses from Ref. the following results are found for 
AtK in the three cases: (I) SAxlQ-'^^sec, (II) S.SxlO-^^sec and (III) GAxlQ-'^^sec. For 
comparison, for neutrino oscillations following pion decay at rest, with Am^ = (leV)^ and 
L = 30 m, Eqn(2.29) gives At^ = 5.6 x 10~^^sec. The result (I) may be compared with 
the mean life of the 6 meson of 1.5x 10~^^sec [1541. Thus the 6 lifetime is a factor of about 



five larger than At^ indicating that Kg — Kl interference should be possible by a similar 
mechanism to neutrino oscillations following pion decays, i.e. without invoking velocity 
smearing of the neutral kaon mass eigenstates. In cases (I) and (II) the interference 
effects observed will depend on the 'characteristic time' of the (non resonant) strong 
interaction process, a quantity that has hitherto not been susceptible to experimental 
investigation 0. If this time is much less than, or comparable to, Atx, essentially equal 
velocities (and therefore appreciable velocity smearing) of the eigenstates will be necessary 
for interference to occur. Since ArriK and Ts are comparable in size, velocity smearing 
effects are expected to be, in any case, much larger than for neutrino oscillations following 
pion decay. These effects are readily calculated using the Gaussian approximation (2.22) of 
the present paper. The main contribution to the velocity smearing is due to the variation 
of the physical mass of the Kg rather than those of the Kl or A. 

For the Bi — B2 oscillation gous to (I) above, T{4S) —* B1B2 {pb = 335 

MeV) the value of AIb is found to be 1.8 x 10~^^sec, to be compared with r(T(4S')) = 



4.7 X 10 sec [0, which is a factor 3.8 smaller. Thus, velocity smearing effects are 



^^A similar physical quantity has been considered in Ref. pql , where the possiblity of observable 
modifications to the exponential decay law and the Breit-Wigner line shape distribution is suggested. 
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expected to play an important role in Bi — B2 oscillations. This is possible, since the 
neutral b-meson decay width (4.3 x 10^^° MeV), and mass difference (3.1 x 10 MeV), 
have similar sizes. 

In closing, it is interesting to mention two types of atomic physics experiments where 



interference effects similar to the conjectured (and perhaps observed |T^, ^7j) neutrino 
oscillations have aleady been clearly seen. 



The first is quantum beat spectroscopy This type of experiment, which has 



previously been discussed in connection with neutrino oscillations |5^, corresponds closely 



to the gedanken experiment used by Heisenberg [15| to exemplify the the fundamental law 
of quantum mechanics Eqn(1.2). The atoms of an atomic beam are excited by passage 
through a thin foil or a laser beam. The quantum phase of an atom with excitation energy 
E* evolves with time according to: exp{—iE*At) (see the discussion after Eqn(3.3) above). 
If decay photons from two nearby states with excitation energies and E"^ are detected 
after a time interval At ( for example by placing a photon detector beside the beam at a 
variable distance d from the excitation foil) a cosine interference term with phase: 

{E:-E*,)d 

(pbeat = = (5.20) 

^atom 



where Vatom is the average velocity of the atoms in the beam, is observed |^8[. An atom 
in the beam, before excitation, corresponds to the neutrino source pion. The excitation 
process corresponds to the decay of the pion. The propagation of the two different excited 
states, alternative histories of the initial atom, correspond to the alternative propagation 
of the two neutrino mass eigenstates. Finally the dexcitation of the atoms and the de- 
tection of a single photon corresponds to the neutrino detection process. The particular 
importance of this experiment for the path amplitude calculations presented in the present 
paper, is that it demonstrates, experimentally, the important contribution to the inter- 
ference phase of the space-time propagators of excited atoms, in direct analogy to the 
similar contributions of unstable pions, muons and nuclei discussed above. 

An even closer analogy to neutrino osillations following pion decay is provided by the 
recently observed process of photodetachement of an electron by laser excitation: the 
'Photodetachment Microscope' A laser photon ejects the electron from, for example. 



an ^^0~ ion in a beam. The photodetached electron is emitted in an S-wave (isotropically) 
and with a fixed initial energy. It then moves in a constant, vertical, electric field that is 
perpendicular to the direction of the ion beam and almost parallel to the laser beam. An 
upward moving electron that is decelerated by the field eventually undergoes 'reflection' 
before being accelerated towards a planar position-sensitive electron detector situated 
below the beam and perpendicular to the electric field direction ( see Fig.l of Ref. |[59|| ). 
In these circumstances, it can be shown that, just two classical electron trajectories 



link the production point to any point in the kinematically allowed region of the detection 
plane. Typical parameters for ^^0~ are [^]: initial electron kinetic energy, 102 keV; 
detector distance, 51.4 cm; average time-of-flight, 117 ns; difference in emission times 
to arrive in spatial-temporal coincidence at the detector plane, 160 ps. An interference 
pattern is generated by the phase difference between the amplitudes corresponding to 
the two allowed trajectories. The phase difference, derived by performing the Feynman 



path integral of the classical action along the classical trajectories |61], gives a very good 
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description of the observed interference pattern. The extremely close analogy between 
this experiment and the neutrino oscillation experiments described in Sections 2 and 3 
above is evident. Notice that the neutrinos, like the electrons in the photodetachement 
experiment, must be emitted at different times, in the alternative paths, for interference 
to be possible. This is the crucial point that was not understood in previous treatments 
of the quantum mechanics of neutrino oscillations. 



Actually, Ref. ||61| contains, in Section IV, a path amplitude calculation for electrons 
in free space that is geometrically identical to the discussion of pion decays in flight 
presented Section 4 above (compare Fig. 3 of the present paper with Fig.3 of Ref. jBTj ) 
The conclusion of Ref. |^ is that, in this case, no interference effects are possible for 
electrons that are mononergetic in the source rest frame. As is shown in Section 4 above, 
if these electrons are replaced either by neutrinos of different masses from pion decay, or 
muons recoiling against such neutrinos, observable interference effects are indeed to be 
expected. 
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Appendix A 

Random thermal motion of the decaying pion in the target has two distinct physical 
effects on the phase of neutrino oscillations, 

<Pl2 (0) = p— + (^1) 

(where the first and second terms in Eqn(Al) give the contributions of the neutrino and 
pion paths respectively): 

(1) The observed neutrino momentum, P^, is no longer equal to Pq, due to the boost 
from the pion rest frame to the laboratory system. (Doppler effect or Lorentz boost) 

(2) The time increment of the pion path to — to (see Eqn(2.16)) no longer corresponds 
to the pion proper time. (Relativistic time dilatation) 



Taking into account (1) and (2) gives, for the neutrino oscillation phase: 
where: 



Pu = 77r-Po(l + cos 61) and 7^ 



Here, ^* is the angle between the neutrino momentum vector and the pion flight direction 
in the pion rest frame. Developing 7^ and v^^ in terms of the small quantity pT^/mT^, 
Eqn(A2) may be written as: 



(j)"^^ (corr) = (1)1^ {0) + 



1 - 



cos 6* 



Ptt 



2Pn 



1 - 



2Pn 



{A3) 



Performing now the average of the interference term over the isotropic distribution in 
cos 9*/. 



(cos 4>i2 {corr))ei 



1 /•! 

-i?e j exp [^012" (corr)] (i cos ^* 



^i?eexp|i0i^''(O) + ( £^ 



1 /Am^L 
6 I Po 





■(• 






1 - 






Po- 



Am^L 

2N 



1 - 



3m^ 
2^ 



(A4) 



In deriving Eqn(A4) the following approximate formula is used: 



1 
2ia 



sm a a 

~ 1 - — 

a D 



where 



a = 



Po 



m, 
P 



< 1 



(^5) 
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The average over the Maxwell-Boltzmann distribution (2.31) is readily performed by 'com- 
pleting the square' in the exponential, with the result: 



-Re exp 



, 2m^Po 



1 _ Hh. 



+ i 



4 \m^J \ Po 



2P^ 



{A6) 



= F^iDop) cos[0^2"(O) + (p^iDop)] 

leading to Eqns(2.32) and (2.33) for the Doppler damping factor F^^Dop) and phase shift 
(t)^{Dop), respectively. 

The correction for the effect of thermal motion in the case of muon oscillations 
(j)'^{Dop), respectively, is performed in a similar way. The oscillation phase: 

{A7) 



2PI 



is modified by the Lorentz boost of the muon momentum and energy, and the relativistic 
time dilatation of the phase increment of the pion path, to: 



mlEf.Am'^L ^ m^Aw?L 



2m,P3 



where 



and Vq is given by Eqn(2.37). Developing, as above, in terms of PTr/m.^, gives: 



(f)i2{corr) 



-"12 



COS 6',,+ — 



2 m,,ArrrL 

Pi 



{A9) 



Performing the averages over ^* and then leads to Eqns(2.49) and (2.50) for the damp- 
ing factor Fi^{Dop) and phase shift (t)^{Dop), respectively. 

The efi'ect of the finite longitudinal dimensions of the target or detector is calculated 
by an appropriate weighting of the interference term according to the value of the distance 
X — Xf — Xi between the decay and detection points (see Fig.l). Writing the interference 
phase as 0i2 = PX, and assuming a uniform distribution of decay points within the target 
of thickness It'- 



(cos 012) = — 



1 

It 
2 



sm ■ 



cos l3XdX 



cos I3L 



/3£t" 2 

= FrargCOSpL (AlO) 

Substituting the value of j3 appropriate to neutrino oscillations yields Eqn(2.34). Since the 
value of 13 is the same for neutrino and muon oscillations, the same formula is also valid 
in the latter case. The same correction factor, with the replacement £t ^ in describes 
the effect of a finite detection region of length in: 



L- 



+ Xi<Xf<L + 



+ Xi 
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Appendix B 



The first step in the derivation of Eqn(4.17) relating Av{fi) to Am^ is to calculate 
the angle d*. in the centre-of-mass (CM) system of the decaying pion, corresponding to 
6 in the laboratory (LAB) system (see Fig. 3). It is assumed, throughout, that the pion 
and muon are ultra-relativistic in the latter system, so that: VT^,Vi{ii) ~ 1. The Lorentz 
transformation relating the CM and LAB systems gives the relation: 

-f^{l + v*{iJ,)cos9t) 

The starred quantities refer to the pion CM system. Making the substitutions: 6*2 = Oi+S, 
0* — 0*j^ S*, Eqns(Bl) may be solved to obtain, up to first order in S, 5* and Am^: 

Av*{f,) = v*M - v*M = 7.(1 + vM cosOirS- vMicosei + v*M)6* (52) 

sin 0^ 

where, (c.f. Eqn(2.37): 

Using Eqn(2.36) Av*{iJ,) may be expressed in terms of the neutrino mass difference: 

[m^ — mj^){m^ + mf^y 
Eliminating now Av*(//) between (B2) and (B4) gives a relation between 5, 5* and Am^: 

^ 7.(l + t;o*(M) cosg*)^^ AmlmlAm'' Ml 

In the LAB system, and in the UR limit, the difference of the velocities of the muons 
recoiling against the two neutrino mass eigenstates is: 

A.(^) = - .,(^) = |M - |M ^ ^[EM - EM] m 

where is the muon energy in the LAB system for vanishing neutrino masses. Making 
the Lorentz transformation of the muon energy from the pion CM to the LAB frames, and 
using Eqns(2.4) and (2.36) to retain only terms linear in Am^ and 5*, enables Eqn(B6) 
to be re-written as: 



Am^ {cos 91 ^v*M) 2 _2^ 



S*{m^ — m^) sin 61 



where E^^ is the energy of the pion beam. By combining the geometrical constraint 
equation for the muon velocities, (4.16) with (B5) and (B7) the angles 5 and 5* may be 
eliminated to yield the equation for LAB frame velocity difference: 

2m^{m^ — m^) B 
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where 



A^{vi{fj.)- v^cos^i) I {cos9l + VQ{fj.)f + 



imlml sin^ 



91 



[m 



m 



^ ^ E^{mi-mf,){l + v*Mcos9i: 



X 



ml{ml-ml) 



imjml sin^ 9^) ^ ^ 

+ (BIO) 



[m: 



To simplify (B8) , the quantity {vi (/x) —v-j^ cos 6*1) is now expressed in terms of kinematic 
quantities in the pion CM system. Within the UR approximation used, 

9i,mT,/E^,m^/E^ < 1 

so that 

- «.oos«. 4 (I - 1| + ^ o((|.)' , (!^)' (mi) 

Writing Eqn(Bl) to first order in 9i, and neglecting terms of 0{9imf): 

' E^{l + v*M^os9l) ^ ' 

Using Eqn(B12), and expressing E^ in terms of pion CM quantities, Eqn(Bll) may be 
written as: 

Expressing the RHS of (B13) in terms of and E^, using the relation: 

, ml(2E^ - E^) - mlE^ 

gives Eqn(4.20) of the text. 

On substituting (B13) into the RHS of (BIO), it can be seen that the factor in the 
large curly brackets is the same in (B9) and (BIO), and so cancels in the ratio A/B va. 
Eqn(B8). It follows that: 

^^(.) - "'l^'"' / cos(?- + r-(/i) \ 

" El{ml-ml) \l + v*Mcos9l) ^""''^ 

Finally, using (B3) and (B14) to express the factor in large brackets in Eqn(B15) in terms 
of E^ and Eqn(4.17) of the text is obtained. 
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Figure 1: The space-time description of — »• i^e oscillations following 7r+ decay at rest, in 
Feynman's formulation of quantum mechanics. In a) a 7r+ comes to rest in the stopping 
target T at time to. The pion, at rest at time to, constitutes the initial state for the path 
amplitudes. In b) and c) are shown two alternative classical histories for the vr"*"; in b),[ c)] 
the pion decays into the mass eigenstate |z/i >, [ |z/2 > ] at times ti, [ ^2 ]• If "^i > ^2, 
and for suitable values of ti and ^2 (^2 > ti), the two classical histories may correspond 
to a common final state, shown in d) where the neutrino interaction e~p occurs 

at time to- As the initial and final states of the two classical histories are the same, the 
corresponding path amplitudes must be added coherently, as in Eqn(1.2), to calculate the 
probability of the whole process. 
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Figure 2: The space-time description of 'muon oscillations' induced by neutrino mass 
differences, following 7r+ decay at rest, in Feynman's formulation of quantum mechanics. 
As in Fig.l, b) and c) show alternative classical histories of the stopped vr^. If rrii > m2 
the velocity of fii is less than that of fi2, and provided that t2 > ti, the muons may arrive 
at the same spatial point at the same time in both classical histories. If the muons are 
detected at this space-time point in any way (not necessarily by the observation of muon 
decay as shown in c)) interference between the correponding path amplitudes occurs, 
according to Eqn(1.2), just as in the case of neutrino detection. 
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Figure 3: Two dimensional spatial geometry for the observation of neutrino or muon 
oscillations following pion decay in flight. Four possible classical histories of a pion, 
originally at the point A, are shown. In the first two, the pion decays either into the mass 
eigenstate |z/i >, at point A or into \i>2 > at point C. If mi > m2, and for suitable values 
of the angles 6i and 62, the neutrinos may arrive at the point B at the same time. If a 
neutrino detection event, such as —>■ e~p, then occurs at B at this time, the amplitides 
corresponding to the paths AB and ACB will be indistinguishable so that they must be 
superposed, as in Eqn(1.2), to calculate the probability of the overall decay-propagation- 
detection process. The third and fourth classical histories are similar, except that the 
neutrino mass eigenstates are replaced by the corresponding recoil muons. The muons in 
the different histories may arrive at point B, at the same time, leading to interference and 
'muon oscillations' if they are detected there. 
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